ApJ in press 

Preprint typeset using I^T^^ style cmulatcapj v. 02/09/03 



THE CMB QUADRUPOLE IN A POLARIZED LIGHT 
Olivier Dore^, Gilbert P. Holder^, & Abraham Loeb^'^ 

^Department of Astrophysical Sciences, Princeton University, Princeton NJ 08544 
^School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 08540 
^Astronomy Department, Harvard University, Cambridge, MA 02138 
ApJ in press 



o 
o 

(N 



(N 



> 
oo 
o\ 

O 

m 
o 

6 

■ 



ABSTRACT 

The low quadrupole of the cosmic microwave background (CMB), measured by COBE and confirmed 
by WMAP, has generated much discussion recently. We point out that the well-known correlation 
between temperature and polarization anisotropies of the CMB further constrains the low multipole 
anisotropy data. This correlation originates from the fact that the low-multipole polarization signal is 
sourced by the CMB quadrupole as seen by free electrons during the relatively recent cosmic history. 
Consequently, the large-angle temperature anisotropy data make restrictive predictions for the large- 
angle polarization anisotropy, which depend primarily on the optical depth for electron scattering 
after cosmological recombination, t. We show that if current cosmological models for the generation 
of large angle anisotropy are correct and the COBE/WMAP data are not significantly contaminated 
by non-CMB signals, then the observed Cf^ amplitude on the largest scales is discrepant at the 
~ 99.8% level with the observed Cf^ for the concordance ACDM model with r = 0.10. Using 
T — 0.17, the preferred WMAP model-independent value, the discrepancy is at the level of 98.5%. 
Subject headings: cosmology: theory - cosmology: observation 



1. INTRODUCTION 

The low quadrupole (and first few multipolcs) of 
the cosmic microwave background (CMB), measured by 
COBE (Bennett et al. 1996) and confirmed by WMAP 
(Bennett et al. 2003), has generated much discussion 
recently, with several papers offering various possible 
causes for suppressed power on very large scales (e.g.. 
Bond 1995, Efstathiou 2003a, Bridle et al. 2003, Con- 
taldi et al. 2003, Tegmark et al. 2003, Chne et al. 2003, 
Feng & Zhang 2003, among others). Often in the case of 
measurements suggestive of new physics, the obvious way 
to advance is to perform a better experiment. However, 
the current measurements of the temperature anisotropy 
power spectrum of the CMB on large angular scales are 
already limited only by how well the Galaxy can be re- 
moved (Bennett et al. 2003), and so the prospects for im- 
proved measurements are poor. In this work we point out 
that polarization measurements on large angular scales 
can test whether the temperature anisotropy on these 
scales is indeed generated within the standard cosmolog- 
ical framework. 

The correlation between temperature and polarization 
anisotropies is well-known (Zaldarriaga & Seljak 1997). 
It has been detected by Leitch et al. (2002) on interme- 
diate angular scales and measured by the WMAP exper- 
iment on the scales of interest here (Kogut et al. 2003). 
This correlation has been used to predict the polarization 
pattern on the sky from the observed temperature pat- 
tern (Jaffe 2003). The correlation between temperature 
and polarization arises because the source of the polariza- 
tion is Thomson scattering of the quadrupole anisotropy 
in the temperature of the radiation field. Spatial fluc- 
tuations in the monopole and dipole of the temperature 
field at the time of recombination seed higher multipole 
anisotropies by free-streaming (for a recent review see Hu 
& Dodelson 2002). Most of the large scale polarized sig- 
nal we consider here originated from Thomson scattering 



of the CMB quadrupole at the relatively recent epoch 
following the reionization of the Universe (Zaldarriaga 
1997). The quadrupole temperature anisotropy seen by 
free electrons at this epoch receives a significant con- 
tribution to its fc-space kernel from density fluctuation 
modes on scales that also contribute to the present- 
day quadrupole (Tegmark & Zaldarriaga 2002). 

For a given realization of the CMB sky, the measured 
temperature anisotropy power at a given multipole, Cf^ 
will have some amount of intrinsic scatter, as will the 
polarization anisotropy power, Cf^ and the cross power 
spectrum, Cf^ . Importantly, these measures of the 
power are correlated. If Cf^ is measured to be low, then 
one would expect a low measure of Cf^ and also a low 
value of Cf^ . Therefore, if measures of Cj^ and Cf^ 
at low multipoles (large angular scales) are not "anoma- 
lously" low, then this would exacerbate the current ten- 
sion with theoretical models (Spergel et al. 2003). 

ACDM is currently the standard cosmological model, 
and so we choose to assume the best-fit WMAP cosmo- 
logical parameters and study how likely it is, from a fre- 
quentist perspective, that the observed data are simply 
realizations of this model. If the assumed model is in fact 
correct, then the correlations between observables pro- 
vide statistical consistency checks. For example, flukes 
of cosmic variance should be partly correlated between 
temperature and polarization observables, so it might be 
expected that large outliers in the temperature data will 
have counterparts in the polarization data. Here, we 
present results of consistency tests using current data, 
and estimate the range of possible future data that could 
be comfortably accommodated within the currently ac- 
cepted cosmological models. 

2. TT-TE CORRELATIONS 

The correlations between the temperature and polar- 
ization power spectra can be used as a powerful consis- 
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Fig. 1. — 3D view (black points) of realizations of WMAP best fit model for £ = 2 — 4, assuming a noise model comparable to one year 
of WMAP observations. Projections are shown on the sides of the "box." Center dark vertical line shows measured WMAP point in this 
space; there is no Cf^ data released at this time so no constraint along the vertical axis. The horizontal ( Cj^ , Cj"^ ) plane corresponds 
to Figure 3 whereas the vertical left ( Cf^ , Cj"^ ) plane corresponds to Figure 4. 



tency test. For example, it is possible to use the mea- 
sured Cf^ to construct a probability distribution for 
Cj^ and compare this to measured values of Cf^ . Al- 
ternatively, the measured Cf^ can be used to estimate 
the most likely values of measured Cf^ or to forecast fu- 
ture Cj^ or Cf^ measurements. We present examples 
of both these calculations below. 

Explicit expressions for the covariance between power 
spectra, assuming Gaussian uncorrelated noise and 
Gaussian beams, are given by Zaldarriaga, Spergel, & 
Seljak (1997) and are useful for understanding the na- 
ture of the covariances. In all that follows we neglect 
beam effects since we are interested in very large scales, 
where £ is much smaller than the inverse of the beam size 
(in radians). All Cf (where X corresponds to TT,TE, 
or EE) are in units of fiK^. 

3. MONTE CARLO METHODOLOGY 

We generate 10^ realizations of the WMAP-only best- 
fit cosmological model aim {^mh^ — 0.13, Q.bh? = 0.023, 
h = 0.68, Us = 0.97, A = 0.8 and r = 0.10 as in table 1 
of Spergel et al. (2003)) and for € < 40 we construct the- 
oretical joint distributions of Cf^ , Cj^ and Cj'^ for 
this model. The scatter in these quantities around the 
input model is both due to cosmic variance and noise, 
and the scatter between different power spectra is some- 
what correlated. We use the 1 year WMAP data for 
Cj'^ and Cj^ from the corresponding publicly-available 
ASCII files 1. 

^ http://laiiibda.gsfc.nasa.gov 



We use a simple noise model which is as close as pos- 
sible to that reported by the WMAP experiment af- 
ter one year of operation. Specifically we assume that 
the white noise level corresponds to the use of 16 W- 
Band, 8 V-band and 8 Q-band channels (see table 1 
in Bennett et al. (2003)) with an effective sky frac- 
tion of 0.86/1.14 (Verde et al. 2003), so that we get 
^ -^ Qy ^ 10-VK per (0.21)^ sq. deg. pixel with 
a corresponding we — wtI^- We add this white noise 
directly to the generated agm- This simple simulation 
scheme neglects the weak extra power at low i originat- 
ing from the residual 1// noise (Hinshaw et al. 2003). 
We neglect any non-gaussian systematic uncertainty that 
may be caused by improper subtraction of the galactic 
component (Bennett et al. 2003; Efstathiou 2003b). 
However, to assess the significance of our assumed level 
of noise, we also consider simulations where the noise 

— 1/2 

amplitude, Wrj^ , was arbitrarily scaled up by a factor 

of \/2, so that the noise contribution to the Cf errors is 
doubled. This excess noise could be either some unknown 
component of low frequency noise or residual galactic 
contamination (although the latter should probably not 
be modeled as a white noise). 

Furthermore, to reduce the effects of covariance due 
to the real cut sky we consider bins of width = 2- 
4, where Cf% are asstmied constant. With such binned 
power spectrum estimates, the covariance between neigh- 
boring bins should be very small and our Monte Carlo 
results should approximate reasonably well the actual co- 
variance of the properly measured Cf (we estimate that 
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Fig. 2. — Realizations of Cj^ and Cj"^ for noise levels typical 
of one year of WMAP observations (2 upper rows) and for cosmic 
variance only (2 lower rows). WMAP data points are marked and 
the measurement noise and cosmic variance are included in the 
realizations (i.e. , it is appropriate for the points to not have error 
bars). The data points in the 2 lower rows are shown only for 
reference (as they contain noise). 



the correlation between our lower bins is less than 10%). 
It is important to note that our procedure is only ap- 
proximate and it is therefore not strictly appropriate to 
compare our results directly to WMAP data; however 
we do not expect that more accurate simulations would 
alter significantly our conclusions. 

4. MONTE CARLO RESULTS 

The result of these realizations is an ensemble of points 
in ( Cj'^ , Cj^ , Cf^ ) space, shown in Figure 1. Slices 
through this volume then provide a frequentist estimate 
of the expected joint distribution of spectra for this par- 
ticular model. In particular, slices passing through ob- 
served points allow investigation of the conditional like- 
lihoods in the other directions. For example, the vertical 
line shows the observed WMAP data and a histogram of 
points along this line would give the conditional proba- 
bility of Cf^ given the observed Cj^ and Cj^ . 

In Figure 2 we show the distribution of expected Cj"^ 
and Cj^ that would be observed, including our approxi- 
mate WMAP noise model. Current data already provide 
a consistency check, even with the relatively large noise 



contribution to Cj^ . On the largest angular scales 
the measurements are approaching the cosmic-variance 
limit for bins of size = 3, even with only one year 
of WMAP data. However, on slightly smaller angular 
scales the data contain a significant noise component, as 
can be seen in Figure 8 of Kogut et al. (2003). This 
noise will be uncorrelated between the Cj^ and Cj"^ 
power spectra and leads to the bins with £ > 8 showing 
little correlation in Figure 2. 

In Figure 2 we also show the Cj^ - Cf^ correlation for 
perfect measurements (no noise) . The current \ow-£ Cf^ 
data are significantly discrepant with the best-fit model, 
but equally remarkable is that the Cf^ measurements on 
the largest scales are not particularly low. The measured 
Cj^ on these scales is close to the middle of the expected 
range if one ignores the correlation with the Cf^ data on 
these scales. However, the probability of the measured 
Cj^ given the observed Cf^ is extremely low if this is 
indeed a realization of the best fit model. The observed 
Cj^ is in fact anomalous by not being low. This can be 
seen comparing the bottom right and top right panels of 
Figure 3. In the bottom right panel it can be seen that 
the observed Cf^ on this scale is just in the range pre- 
dicted by the best-fit model (26% of the models lie above 
the observed value). In the top right panel we see that 
this agreement disappears when we apply the condition 
that the observed Cf^ is low (0.17% of the models lie 
above the observed value (or 0.85% when the noise 
is scaled up arbitrarily by a factor of 2). In realizations 
of the best-fit model it is rare that Cf^ is as low as 
the value measured by WMAP, but in the few realiza- 
tions where Cf^ was low it was usually the case that 
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Fig. 3. — Lowest multipolc bin distribution of Cj"^ and 
assuming the WMAP best fit cosmological model and one year of 
observation (top left). The bottom right panel shows the likelihood 
distribution of Cf^ with no information on Cf"^ , and the top right 
panel shows the likelihood subject to the constraint that Cf^ has 
the measured value. The bottom left panel shows the likelihood of 
the Cj"^ given the observed Cf^ . Dashed vertical lines indicate 
the measured WMAP values. The dashed distributions correspond 
to simulations where the noise level has been arbitrarily increased 
by a factor \/2- 
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Fig. 4. — Same as Figure 3, now for Cf^ - Cf^ covariance, again 
assuming a noise model comparable to one year of WMAP data. 
Here we show the multipole bins 2 < £ < 4 (left) and 20 < £ < 23 
(right), the two largest outliers in the Cj"^ bins considered. 



the Cj^ was also low. Correspondingly, middle-of-the- 
road Cj^ values are unlikely to appear with low Cf^ 
, as shown in the lower left panel. Specifically, 99.9% of 
the Cf^ values are greater than the measured one given 
the measured value of Cf^ , while when the measured 
Cf^ is not included the fraction drops to 99.0% (with the 
corresponding numbers being 99.8% and 99.1% when the 
Cf noise is arbitrarily doubled). These results are dic- 
tated by the profile of the joint ( Cf^ , Cj^ ) distribution 
in the top left panel of Figure 3. 

Upcoming measurements of the Cf^ power spectrum 
on large scales could shed some light on this problem. 
At the noise levels expected in the near future there is 
little correlation between the Cf^ and Cf^ power spec- 
tra, as shown in Figure 4, where the histograms in the 
right top and bottom panels are nearly unaffected when 
the temperature information is included. The 3D plot 
in Figure 1 again shows that the observed Cj^ - Cf^ 
large angle pair (shown as a vertical line) is exceedingly 
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Fig. 5. — Probability for Cj'^ larger than the observed value, 
given the observed Cf^ value for 2 < Z < 4, V{Cj^ > Cjf^^\Cj^), 
as a function of the optical depth value r (solid line). Also shown, 
for reference, is the same probability without the condition on the 
observed Cj'^ value (dashed line). All cosmological parameters 
arc kept the same as in our fiducial ACDM model except for t and 
ris that are changed simultaneously along the t-Us degeneracy line 
favored by WMAP. The power-spectrum amplitude is marginalized 
over. Note that P(Cj-^ > Cjg'^\CfJ') does not exceed 5%. 



unlikely, but also shows that there is a fairly tight cor- 
relation between Cf-^ and Cf^ on these scales. This 
provides an important consistency check on our under- 



standing of the CMB. The observed C'[^ - Cj"^ pair on 
these scales reduces significantly the expected range of 
Cf^ (assuming the best fit model). 

5. COSMOLOGICAL IMPLICATIONS 

The correlations between power spectra have several 
cosmological implications. 

So far we assumed the best-fit WMAP value for the 
optical depth to electron scattering after cosmological 
recombination, t = 0.10 (Kogut et al. 2003; Spergel et 
al. 2003) and a power-law primordial power spectrum. 

Adopting a different optical depth r has a signifi- 
cant but still limited effect on our results. We repeated 
the analysis by considering various optical depth val- 
ues and the results are shown in Figure 5. The solid 
line shows the probability for 2< i <4 of having Cj^ 
larger than the observed value given the observed Cj"^ 
value, V{Cj^ > Cj^^\Cj^), and the dashed hue shows 
the same without the condition on the observed Cj'^ 
value, V{Cj^ > Cj^). All cosmological parameters were 
held fixed except for t and ris which were varied in the 
ranges 0,05-0.29 and 0.95-1.05 respectively, along the 
Us — T degeneracy line favored by WMAP [see Figure 5 
of Spergel et al. (2003)], The power spectrum amphtude 
As was marginalized over. Other parameters will have 
little impact on the large-angle polarization, as shown in 
Kaplinghat et al. (2003), A complete frequentist treat- 
ment, allowing all parameters to vary from their best fit, 
is unlikely to lead to qualitatively different results. 

Increasing r makes the observed Cj^ more likely, 
but is not sufhcient to alleviate the tension between the 
low Cf^ and the average Cf^ : note that V{Cj^ > 
Cjfp\Cj^) does not exceed 5%, Note that this tension is 
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Fig. 6. — Same as Figure 3, except assuming cosmic variance 
limited measurements and a complete suppression of the primordial 
power spectrum at wavenumbers fc < fc* = 3x lO"** Mpc"'^ (Cline, 
Crotty, & Lesgourgues 2003). 



not included in current estimates of the optical depth. In 
previous analyses of WMAP data, the Cf^ likelihood is 
multiplied by the bottom right panels rather than the top 
right panels of Figure 3. For all but the largest angular 
scales (the first few multipoles) this effect is negligible, 
but it is clear from Figure 5 that the direction of the bias 
is such that the current estimates of r are likely to be 
low. 

Additional suppression of the matter power spectrum 
on large scales tends to reduce the discrepancy of the 
current Cf^ data, but not by a large amount. The 
mapping from matter power spectrum to temperature 
anisotropics is fairly broad in Fourier space (Tegmark & 
Zaldarriaga 2002), so that suppression of power on the 
largest scales in the matter power spectrum does not lead 
to a sharp suppression only on the largest angular scales 
in the CMB, which is what the data seems to suggest. 
At the same time, the Cf^ multipoles at £ < 5 will 
be suppressed as well (Bridle et al. 2003; Chne, Crotty, 
& Lesgourgues 2003), making it even more difficult to 
match the observed Cj^ - Cf^ pair as illustrated in 
Figure 6. Whereas now only 3.2% of the models lie above 
the observed Cf^ value for our lowest £ bin (as compared 
to the previous 26%), less than 0.02% (as compared to 
0.17%) lie there if we include the measured Cf^ . 

6. CONCLUSIONS 

The correlations between the temperature and polar- 
ization power spectra provide a powerful consistency 
check on CMB anisotropy measurements. Cosmic vari- 
ance fluctuations should be correlated between power 

^ http:/ /www. astro. esa.int/SA-general/Projects/Planck/ 



spectra. By assuming the best-fit cosmological model 
we have shown that the observed Cf^ - Cf^ data at 
2 < £ < 4: are alarmingly large outliers. Much has been 
made recently of the degree to which the Cf^ data on 
these scales is anomalously low, but it is almost equally 
alarming that the Cf^ is low and Cf^ is apparently typ- 
ical. In most realizations of the WMAP best-fit model 
with low Cf^ on large scales, Cf^ is also low. It is 
extremely unlikely to high as the measured 

value, given the low observed Cj"^ . This correlation 
is not currently included in likelihood analyses of CMB 
data, but should be relatively easy to incorporate in fu- 
ture work. Current estimates of the optical depth are 
likely biased low. Prescriptions for reducing the primor- 
dial quadrupole may have problems producing Cj^ am- 
plitudes as high as those that are observed on the largest 
angular scales, given the correlation between Cf^ and 
rTE 

Most of the other multipole bins up to ^ < 20 in Figure 
2 appear reasonably consistent with the best fit cosmo- 
logical model. The £ = 20 — 23 bin is discrepant at nearly 
the same level in the Cj^ power spectrum, but at cur- 
rent noise levels the degree of correlation between the 
measured Cj^ and Cj^ is negligible. The full four-year 
data of WMAP will help to improve these constraints. 
More accurate treatments of the errors and residual cor- 
relation due to the cut-sky will be possible. Furthermore, 
it would be particularly interesting to study the value of 
the probability V{Cj^\Cj^) for various galactic cuts as 
a probe of the galactic contribution (D. Spergel, private 
communication). In the more distant future, data of the 
quality forecasted for the Planck satellite^ will provide 
even more powerful consistency checks on the best-fit cos- 
mological model. More accurate modeling of the galactic 
emission will also certainly help to address those issues. 

What would it mean if the observed statistics of the 
CMB do not appear to be consistent with the best fit cos- 
mological model? The simplest explanation would be 
that at least one of the measured components is not 
purely cosmological in origin, possibly due to Galactic 
contamination. For example, removing the Galactic fore- 
ground appears to enhance the inferred quadrupole (de 
Oliveira-Costa et a/. 2004). The absence of foreground 
contamination would be an exciting indication of new 
physics. 
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